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Abstract. We consider a two-mode atomic Josephson junction real- 
ized with dilute dipolar bosons confined by a double-well. We employ 
the two-site extended Bose-Hubbard Hamiltonian and characterize the 
ground-state of this system by the Fisher information, coherence visibil- 
ity, and entanglement entropy. These quantities are studied as functions 
of the interaction between bosons in different wells. The emergence of 
Schrodinger-cat like state with a loss of coherence is also commented. 



1 Introduction 

Ultracold and dilute bosons confined in one-dimensional (ID) double-well poten- 
tials [U are the ideal framework to address the study the Josephson effect [5], and 
the formation of macroscopic coherent states |3I4I5I6I7] and Schrodinger-cat states 
[819110 11 12 14 15 . The coherent dynamics of the bosonic cloud in the double-well 
potential (bosonic Josephson junction) |1I16| . widely studied with alkali- metal atoms, 
is very well described by Josephson equations [3], and their extensions |17ll8ll9l^ . 
The Josephson equations are valid when the inter-atomic interaction is weak against 
the Josephson coupling energy (i.e., the tunneling probability amplitude multiplied 
by the number of bosons) and the number of bosons is much larger than one. In this 
case, semiclassical approximations can be performed and the system is considered in a 
coherent state [3T] . By increasing the coupling strength of the boson-boson repulsive 
interaction |4I5I6I7| the crossover from a coherent state (superfluid-like regime) to a 
pure Fock state (Mott-like regime) takes place. For attractive bosons, the Josephson 
equations predict the spontaneous symmetry-breaking when the strength is above a 
critical value [3120] . while the two-site Bose-Hubbard model [22] predicts the forma- 
tion of a Schrodinger-cat state I8I9I10I11I12I14I15] . When the attraction between the 
bosons becomes strong enough, one expects the collapse of the cloud 2 23) . 

The above considerations hold for bosons with negligible dipole or electric mo- 
ments, so that the dipole-dipole interatomic interaction can be be safely neglected. 
This is not the case for junctions made by dipolar bosons, as for example ^^Cr, char- 
acterized by very large atomic magnetic dipoles. 

In the present contribution we study the emergence of cat-like states in atomic 
Josephson junctions made of dipolar bosons. Recently, dipolar quantum gases, where 
considering long-range and anisotropic dipole-dipole interaction between magnetic 
or electric dipoles makes sense, have attracted a lot of interests ^24,25) . To date, 
important theoretical efforts have been devoted in order to investigate dipolar bosonic 
gases trapped in double |26l27l28j and triple wells [29]. Remarkably, Abad and co- 
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workers [35] have shown that the dipolar interaction makes possible self-inducing a 
double-well potential structure. 

We use as theoretical tool the extended two-site Bose-Hubbard (EBH) Hamilto- 
nian .29 30 31] where both the on-site (intra- well) interaction and nearest-neighbor 
(inter- well) density-density one are considered. We follow the same path as in [13114115] 
Thus, we diagonalize the two-site EBH Hamiltonian and study the Fisher informa- 
tion F, the coherence visibility a, and the entanglement entropy S of the ground-state 
by fixing the on-site interaction amplitude and varying the inter-well density-density 
interaction. We find that the presence of a macroscopic superposition state corre- 
sponds to a sufficiently large values of the Fisher information F, as expected from 
the studies presented in |33l34j . By increasing the inter- well interaction, F grows, the 
coherence visibility goes to zero, and the entanglement entropy S reaches a maximum 
value - in correspondence to a given inter- well interaction. For larger nearest-neighbor 
interaction, S decreases and the ground-state of the EBH evolves toward a cat-like 
state. 

As a second task, we solved the ordinary differential equations for a bosonic junc- 
tion obtained by using the quasi-classical coherent state, varying the inter-well inter- 
action. When such an interaction becomes strong enough, the population imbalance 
starts to oscillate around a non-zero value. The on-set of this behavior takes place for 
that inter- well interaction corresponding to the maximum value of the entanglement 
entropy S. 



2 The model Hamiltonian 

We consider an ultracold gas of N identical dipolar bosons of mass m. The boson- 
boson interaction derives from the sum of a short-range contact potential and a 

2 cos*^ 

long-range dipole-dipole potential, i.e. V{r — r') = gd{r — r') + 7 — ■ — , where 

|r — r'l'' 

g = Anh as/m with Og the interatomic s-wave scattering length; 7 = ^Qfj? /An for 
magnetic dipoles {fio is the vacuum magnetic susceptibility and n is the magnetic 
dipole moment) or 7 = (P/Aweo for electric dipoles {eq is the vacuum dielectric con- 
stant and d is the electric dipole moment). We assume that the bosons are polarized 
by a sufficiently large external field with all the dipoles aligned along the same direc- 
tion; 9 is the angle between the vector r — r' and the dipole orientation. We suppose 
the bosons are confined by the superposition of an isotropic harmonic confinement 
in the transverse (y — z) radial plane and a symmetric double-well potential Vdw in 
the axial direction (x): Vtrap(r) = Vowix) + ^muj'^{y'^ + z^) with wj, the trapping 
frequency in the radial plane. In the following we shall consider a strong transverse 
confinement so that the system can be treated as one-dimensional (ID). In particular, 
the transverse energy Tiu)^_ is assumed much larger than the characteristic energy of 
bosons in the axial direction. Then, a dilute bosonic gas can be described with the 
help of the two-mode Bose-Hubbard Hamiltonian [52]. We are considering particles 
that interact with each other via long-range forces; thus the interactions between 
bosons in different wells have to be taken into account. Proceeding from the second 
quantized Hamiltonian written in terms of the space dependent bosonic field opera- 
tors, one integrates out the spatial degrees of freedom, namely r and r', involved in 
the trapping and boson-boson interaction potentials and in the kinetic term as well. 
Thus one obtains the following two-site version of extended Bose-Hubbard (EBH) 
Hamiltonian [3T] (in the following L will stand for left, R for right): 

H = - j{a\aB. + o^b^cll) + ^{fiLinL - 1) + nninR - 1)) + UxfiLriR , (1) 
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where afe,a|, (fc — L,R) are bosonic operators and fik = a^dk counts the number of 
particles in the fcth well. Notice that due to the above mentioned integration over the 
spatial degrees of freedom, the quantities J, Uq, and Ui do not depend on the spatial 
coordinates, but only on the microscopic parameters of the system, see, for example 
[19] . In particular, J is the hopping amplitude between the two wells, [/q the on-site 
interaction amplitude, and Ui the nearest-neighbor interaction amplitude. 

3 Analysis 

In the absence of the nearest- neighbor interaction {Ui — 0), the ground-state of 
the Hamiltonian ([1]) in the limit Uq/ J ^ —oo is the macroscopic superposition state 
-^{\N,0) + \0,N)), also known as "NOON" state or "Schrodinger cat state" . However, 

when the interaction between the bosons is attractive, i.e. < (resulting mUo/J < 
0), the collapse of the atomic cloud can be observed if the atomic density is sufficiently 
high. Apart the question of the possible collapse, the realization of the cat state is not 
trivial due to the very tiny separation (in the presence of finite couplings) between 
the two lowest levels that makes the cat state very fragile, see, for example [TU] . 

Let us suppose that Ui ^ and J = 0. It is easy to see that, for Ui > Uq, the 
Hamiltonian ([T]) admits the NOON state as its lowest eigenstate even if the on-site 
interaction is repulsive. This allows to realize cat-like states by-passing the collapse 
problem. It is worth to observe that, in the absence of the hopping, the emergence of 
the NOON state can be understood in terms of the interplay between the intra-well 
interaction, in practice Uofil {k = L,R), and the inter-well one, i.e. UifiLfiR. In fact, 
the effect of the former is to establish a balanced population among the two sites, 
while the latter is minimum when one of the on-site average occupations vanishes. 

We therefore study the emergence of cat-like state by fixing Uq and varying Ui for 
a given N. This is possible since - see, f.i., - the on-site interaction results from 
short-range interaction and dipole-dipole one, while Ui depends only on the dipole- 
dipole interaction: then, we have to vary both contact and dipole-dipole interaction 
simultaneously, in such a way that we can change Ui keeping Uq constant. 

Since the total number of bosons is a preserved quantity, we numerically solve the 
eigenproblem H\Ej) = Ej\Ej) for a fixed number N of bosons. For each eigenvalue 
Ej {Eq < El < E2 < ...) with j = 0, 1, N, the corresponding eigenstate \Ej) will 
be of the form 

N 

\E,)=Y,c^^\^,N-^) . (2) 

1=0 

Let us fix the rescaled on-site interaction U = Uq/J ai a positive value and vary the 
rescaled nearest-neighbor interaction V = Ui/J. We find that it exists a crossover 
value of V, say Vcr- When V exceeds such a value, the system begins to lose the 
coherence. This crossover connects a ground-state with the maximal probability at 
i = N/2 (single peak centered around i = 0) - for = - to a ground-state which has 
the maximal probability at i = n and i = N — n (two separated peaks symmetric with 
respect to i — 0) with n approaching as F increases (this last situation corresponds 
to the emergence of the Schrodinger cat state) , as it can be clearly observed from the 

panel a) of Fig. 1, where we have plotted the coefficients |c|°^p. We have numerically 
found that the V crossover value is related to U by the following linear relation: 

Vcr = U + A, (3) 

where A > depends on N. For example, when N = 30, A — 0.075, while when 
N = 100, A is equal to 0.0212. Translated into the language of the atomic Josephson 
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Fig. 1. (Color online), a) Coefficients |c^''''[^ for the ground-state as a function of i/N for 
U = 0.2. b) The normalized quantum Fisher information F vs V ai U = 0.2 (solid line) and 
U — 0.8 (dotted line), c) The coherence visibility avsV at U — 0.2 (solid line) and U = 0.8 
(dotted line), d) The entanglement entropy S ys V at U — 0.2 (solid line) and U = 0.8 
(dotted line), e) and f) The population fractional imbalance z(t) vs time t at U — 0.2. Initial 
conditions 2(0) = 0.1, 0(0) — 0. Ic'"']^, i/N, U, and V are all dimensionless quantities. 
Number of bosons A*' = 100. Time is in units of h/J. 



junctions: when V approaches Vcr the junction is entering the self-trapping regime in 
which the fractional population imbalance between the two wells, z = {Nl — Nfi)/N 
{Nk is the number of bosons in the fcth well) oscillates around a non-zero time averaged 
value, see the panel f) of Fig. 1 and the discussion below. 

From an analytical point of view, the above quoted interpretation can be supported 
by calculating the expectation value E of the EBH Hamiltonian ^ with respect to 
the state [H] \QC) = \CS)l'S> \CS)r, where \CS)k {k = L,R) - which describes the 
Bose-Einstein condensate (EEC) in the fcth well - is [ID] 

(4) 

The complex quantity Zk is the eigenvalue of the annihilator in the fcth well, i.e. 

ak\CS)k^Zk\CS)k. (5) 
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The absolute values of the two Zk are related to the average occupancy of the two 
wells: 

Nk = {QC\hk\QC) = Izfcp. (6) 

Therefore Zk are conveniently parametrized as z^ = V-^fe exp(i0fc), where 6k are phase 
variables. Following the same procedure as in when only Uq was present, we get 
for the expectation value E the following expression: 



S= (QC|7?|QC) = -TV JVT^cos 0+^^(1 + + ^ (7) 

where (j) — On — Ol {Ok is the phase of the BEC in the fcth well) and the fractional im- 
balance z = {\zl\^ — \zr\'^)/{\zl\'^ + \zji\'^). Starting from the energy ([7]) , the equations 
of the motion for the generalized coordinate (f> and its conjugate moment z provide 
the ordinary differential equations (ODEs) for z{t) and (pit). These ODEs, with the 
time scaled with respect to fi/J, read 



z{t) = - Vl - z^{t) sin</>(i) 

U — V^^,. z(t) cosdiit) 
Ht)^^^Nzit)+ ) (8) 

By solving these ODEs, we have obtained the panels e) and f) of Fig. 1. Now, it is 
straightforward to show that the non-zero minimum of the energy ([T]) is obtained 
with 6 = and 



It is easy to show that stationary points of the form (z,0) (with z 7^ and zero 
relative phase 0) for the ordinary differential equations (O can exist only when V > U 
(Ui > Uq). In particular, the point (zsb,0), see Eq.®, is a stationary point for the 
afore mentioned ODEs if the order relation in Eq.® is met. Notice that if the ODEs 
^ are solved with the initial condition (j){t = 0) = 0(0) = 0, oscillations of the 
population imbalance around a non-zero time averaged value {{z(t)) 7^ 0) are possible 
for any z{t = 0) = z(0) provided that V > U + 2/N [Ui > Uq + 2J/N). For this 
kind of oscillations, the temporal average of the relative phase is zero. On the other 
hand, if C/ > V - i.e. Uq > Ui - (and hence, in particular, in the ordinary BJJ case 
corresponding to Ui =0), for 0(0) = 0, a non-zero z cannot be a stationary point 
for the above ODEs. In this case, the oscillations around {z{t)) ^ take place if z(0) 
exceeds a critical value and the phase is an increasing function of the time, and we 
have the so called running phase mode [3]. 

At this point it is worth to comment about the following. The z-symmetry bro- 
ken states that one observes around stationary stables point of ([5]) are related to the 
semiclassical limit that we have used above. Such states have a large lifetime that 
scales exponentially with A'^ [3]. In a full quantum two- mode approximation, instead, 
the ground state is always symmetric in the population imbalance, being well approx- 
imated by a symmetric superposition of two unbalanced coherent like states, namely 
a Schrodinger cat state 15 . 

We see that at varying U (that is, the intra- well interaction), V = U + 2/N - that 
is nothing but Eq. ^ - gives the value of V (that is, the inter-well interaction) which 
signals the onset of the self-trapping regime in the bosonic Josephson junction. 

By following [14 , we characterize the ground-state of our system by studying three 
indicators: the quantum Fisher information, the coherence visibility, and the entangle- 
ment entropy. We shall analyze these quantities as functions of V in correspondence 
to the same on-site rescaled interaction. 
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It is useful to recall the definition of the quantum Fisher information (QFI) that 

we shall use in the following: Fqpj = {(jil — fiR)'^) — {{^l — "-i?)) |37j . The quantum 
Fisher information is one of the key quantities of the estimation theory |38| . In par- 
ticular, the QFI is related to the bound on the precision with which the accumulated 
phase in an interferometric experiment can be determined. For a given input state 
of an interferometric procedure, the best achievable precision is the reciprocal of the 
squared root of QFI, known as the quantum Cramer-Rao bound [35] • Notice that the 
definition of QFI above provided is relevant for an atomic interferometer based on ro- 
tations of the input state about the z axes of the Bloch sphere. However, other choices 
would be possible, namely those relying on input state rotations about the x and y 
axes of the Bloch sphere [34137) . Moreover, the quantum Fisher information is a pa- 
rameter related to the multiparticle entanglement, namely to the indistinguishability 
of the bosons [31] . 

It is convenient to normalize Fqpj at its maximum value N"^ - as well known from 
[33134] - by defining the normalized quantum Fisher information F = Fqfi/N^- In 
terms of the coefficients cf^ , F is given by: 

^-]^E[2^-^]'K^°^I'- (10) 

In the large N limit we can compare our results to analytical asymptotic behaviors 
[T5] . For V > U + 2/N, in fact, we expect to have 

F{V)c^^{l + iN-l)zl) , (11) 

with Zsb given by Eq. ((H]). In the extreme NOON state {zsb = 1) we obtain F = 1. 
For U > V, instead, we get the following large N limit result 



Let us move, now, to the coherence visibility a given by a = 2 \{a[aR)\/N [3] which 
characterizes the degree of coherence between the two wells. The expectation value of 
the operator a'^aR is evaluated in the ground-state \Eo) and the visibility a is given 

by 

"=^E-rM?iV(* + l)(A^-0. (13) 

i=0 

In the large N limit and for V > U + 2/N, the asymptotic behavior is given by 

aiV) ^l~zl , (14) 
while its value for U > V , always for large number of bosons, is well approximated 

by m 

Finally, we calculate the entanglement entropy S [39. . This is an excellent measure 
of the quantum entanglement of the ground-state \Eq). When the system is in \Eo), 
the density matrix is p = |_Eo)(£'o|- S is defined as the von Neumann entropy of the 
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reduced density matrix pl — Trjip, which is the matrix obtained by partial tracing 
the total density matrix p over the degrees of freedom of the right well. S measures the 
bi-partite entanglement, that is the amount of genuine quantum correlations between 
the left well and the right one, seen as two partitions of the whole quantum system. 
The entanglement entropy S is given by 

N 

S = -TrpL log, PL^-J2 l^r' P log2 Icf ^ r ■ (16) 

1=0 

We observe that the cat state is not the maximally entangled state - from the left- 
right bi-partition perspective - for our system. In fact, from Eq. pB]) we see that, for a 
NOON state, S* = 1, while the maximum value attained by the entanglement entropy 
for = 100 is Smax — 5.4, see panel d) Fig.[TJ From analytical considerations based 
on the asymptotic shape of the reduced density matrix [15l41j we find that the peak 
of the entropy weakly depends on the interacion parameters and, in the large N limit, 
is simply given by 

Smax ^ ^\og2i2neN) , (17) 

while the location of the peak is around V = U + 2/N, in perfect agreement with the 
numerical results. The value of at F = 0, in the large N limit, is recovered [T^ 

^(^ = 0)4log,(.e/^^), (18) 

as well as its asymptotic value ior V oo which is S' — )■ 1, sign of the NOON state. 

f , a, and S', as functions of V are shown in Fig. 1. We see that in correspondence 
to an increasing of the nearest-neighbor interaction the ground-state of the system 
experiences a "catness" enhancement accompanied by a softening of the coherence 
visibility and of the entanglement entropy which gets its maximum value when V = 
Vcr- This trend can be retrieved also for larger values of U. For each of three indicators 
- F, a, S - in Fig. 1, we show the comparison for two different values of U. We have 
found - as expected from Eq. ^ and the following discussion - that an increasing of the 
repulsive interatomic interaction produces a shift of Vcr towards larger values. Notice 
that also when the cross-well interaction is absent, = 0, the entanglement entropy 
S attains its maximum value when the ground-state of the underlying Hamiltonian 
is a cat-like state |14ll5l42j . 

By means of long-range potential one could circumvent the collapse - and thus 
the instability - of the bosonic cloud which should take place with sufficiently high 
densities for attractive interactions. Indeed, the long-range potential between the 
confined bosons introduces an additional degree of freedom with respect to the short- 
range one - that is the nearest-neighbor density-density interaction - which when 
suitably tuned makes possible to get the NOON state even for repulsive interactions. 
We stress that such a phenomenon is absent in the standard two-site Bose-Hubbard 
model, i.e. in the absence of inter- well interactions [10113114115142143] . 

4 Conclusions 

We have considered an atomic Josephson junction made of a finite number of interact- 
ing dipolar bosons. By employing the extended two-site Bose-Hubbard Hamiltonian, 
we have carried out the zero-temperature analysis by finding the ground-state of the 
system and characterizing it with the help of three indicators: the Fisher information. 
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the coherence visibility, and the entanglement entropy. We have studied these quan- 
tities by fixing the on-site interaction and varying the nearest-neighbor one. We have 
found that the presence of a Schrodinger-cat like state in the double- well corresponds 
to sufficiently large values of the Fisher information. We have pointed out that the 
cat-like state emerges when, within the underlying classical model for the junction, 
the population imbalance parameter oscillates around a non-zero value. This kind of 
oscillations sets in for the nearest-neighbor interaction signing the maximum of the 
entanglement entropy. In this situation the coherence visibility is quite small but it 
increases as the inter-well interaction strength becomes sufficiently small. 

As for what concerns the future perspectives, the bosonic dipolar interaction might 
be very useful in dynamical generating cat-like states with one L44_ and two bosonic 
components |45j since it is able to make wider the separation with the first excited 
state, making in this way the cat state more robust against decoherence. 

The present work has been supported by Progetto Giovani (University of Padova): 
" Many Body Quantum Physics and Quantum Control with Ultracold Atomic Gases" 
and by Progetto di Ateneo (University of Padova) : " Quantum Information with Ul- 
tracold Atoms in Optical Lattices" . 
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